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It has recently been predicted that certain magnetic materials contain mobile 
magnetic charges or 'monopoles' 1 . Here we address the question of whether these 
magnetic charges and their associated currents ('magnetricity') can be directly 
measured in experiment, without recourse to any material-specific theory. By 
mapping the problem onto Onsager's theory of weak electrolytes, we show that 
this is possible, and devise an appropriate method. Then, using muon spin 
rotation as a convenient local probe, we apply the method to a real material: 
the spin ice Dy2Ti2C>7. Our experimental measurements prove that magnetic 
charges exist in this material, interact via a Coulomb interaction, and have 
measurable currents. We further characterise deviations from Ohm's Law, and 
determine the elementary unit of magnetic charge to be 5 /xbA \ which is equal 
to that predicted by Castelnovo, Moessner and Sondhi 1 using the microscopic 
theory of spin ice. Our demonstration of magnetic charge transport has both 
conceptual and technological implications. 



1 Introduction 



The transport of electrically charged quasiparticles (based on electrons, holes or ions) plays 
a pivotal role in modern technology as well as determining the essential function of bio- 
logical organisms. In contrast, the transport of magnetic charges has barely been explored 
experimentally, mainly because magnetic charges are generally considered to be, at most, 
convenient macroscopic parameters^ rather than sharply defined quasiparticles. However, 
the recent proposition of emergent magnetic monopoles in spin ice materials^ may change 
this point of view. It is of great interest to investigate the transport of such atomic-scale 
magnetic charges, as the resulting 'magnetricity' could have potential technological signif- 
icance. Furthermore, the direct observation of magnetic charge and current would afford 
the strongest case for believing in their reality. Thus, in the analogous case of electrical 
charge transport (electricity), the fact that both current / and the elementary charge e can 
be directly measured emphasises that these quantities are 'real', rather than just convenient 
parameters in a model. 

Spin ices are frustrated magnets, a class of magnetic material that is well known for sup- 
porting exotic excitations, both in experiment^ and in theory^. The spin ice family^^ 
is predicted to support sharply defined magnetic monopole excitations and offers the pos- 
sibility of exploring the general properties of magnetic charge transport in an ideal model 
system^. There is significant experimental evidence to support the existence of spin ice 
magnetic monopoles^^, but this evidence does not include the direct observation of charge 
or current and relies strongly on interpretation via the microscopic theory of spin ice. In the 
analogous electrical case, charge and current may be directly measured, with only the basic 
knowledge that the sample is say a semiconductor or a metal. It therefore seems realistic to 
seek a method of measuring magnetic charge and current that is similarly direct and robust. 

As spin ices are magnetic analogues of water ice, the magnetic charges of spin ice are 
predicted to be directly analogous to water ice's mobile ionic defects. Water ice in turn may 
be classified as a weakly dissociated electrolyte, for which there exist a variety of general 
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experimental methods to prove the existence of unbound ions and to estimate their charge 
and mobility. Although most of the relevant methods have no possible magnetic analogy, 
we have identified one that does. The 'second Wien effect' describes the nonlinear increase 
in dissociation constant K (or equivalently the conductance) of a weak electrolyte (solid or 
liquid) in an applied electric field $ n l 12 [ In a seminal work of 1934 Onsager^ derived a 
general equation for the second Wien effect. This equation provides an excellent description 
of experimental conductivity measurements, and remarkably for a thermodynamic relation, 
enables the determination of the elementary charge e. 

The aim of the current Article is to present a theory for the direct characteristion of 
magnetic currents via the Wien effect, and to test it experimentally. The principal theoretical 
predictions are that the magnetic conductivity is proportional to the fluctuation rate v of 
the magnetic moment, and that the 'elementary' magnetic charge Q (which should depend 
on material) may be derived from the initial slope and intercept of the field dependence of 
v(B) via the equation: 

Q = 2.1223 m 1/3 T 2/3 . (1) 

Here m = slope/intercept is the field gradient of the relative magnetic conductivity, T is the 
temperature and the tilde means that Q is measured in units of /U#A 1 (SI units are used 
elsewhere). Eqn. [l] is valid for a sufficiently weak magnetic electrolyte and assumes a small 
relative permeability of the material (although this assumption is not strictly necessary: see 
below). 

Thus we predict a specific data collapse of v(B, T) over the range of field and temperature 
where the theory is valid. As the charge depends only on the ratio u(B)/u(0), we simply 
require an experimental quantity that is proportional to v{B) in the weak field limit. This use 
of relative quantities obviates the need for absolute measurements, which makes the method 
particularly robust and flexible. In the case of the spin ice Dy2Ti 2 07'^', we have found that 
transverse field muon spin rotation^ affords a convenient probe. We have therefore used 
this technique to test the theory and to determine the elementary magnetic charge in the 
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spin ice Dy 2 Ti 2 07. 

The Article is organised as follows. In the "Theory of the Magnetic Wien Effect" Section, 
we describe the magnetic equivalent of the second Wien effect, and derive Eqn. [T] In the 
"Application to a Real Material" Section we describe our muon spin rotation (/zSR) exper- 
iments on Dy 2 Ti 2 07. Finally, in the "Discussion and Conclusions" Section, we summarise 
our main findings and discuss the broader implications of the work. 



2 Theory of the Magnetic Wien Effect 

The dissociation of a weak electrolyte may generally be represented by two successive equi- 
libria, the first representing the formation of an associated or bound ion pair and the second 
the dissociation into free ions. For example, in the specific case of water ice, the equilibria 
may be represented: 

2H 2 = [HgO+OIT] = H 3 + + Otr, (2) 

where [. . .] represents the bound pair. Here the unbound ions interact according to a (pos- 
sibly screened) Coulomb law. A physical picture of the Wien effect is that an applied field 
accelerates the free ions and, opposed by Brownian motion, in some cases does enough work 
to overcome the Coulomb potential barrier that binds the ions together. The result is an 
increase, with field, of the rate of dissociation and hence of the corresponding equilibrium 
constant. The field acts only on the forward reaction of the second (dissociation) equilibrium 
via the electrical force F = —eE. Onsager's theory^ is valid under the condition that the 
concentration of unbound defects is sufficiently small for the Debye screening length to be 
much greater than the association distance 1 12 1 14 1 . His main result may be approximated by 
the following form in the weak field limit: 

K{E) = K{0)(l + b + £..), (3) 
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valid for b < 3, where b is the dimensionless group 

e 3 E 



8nek 2 T 2 



(4) 



Here the symbols have their usual meaning and e is the permittivity of the solvent. The 
quantity —kTVSb may be interpreted as the Coulombic barrier to ion pair dissociation, at 
which the field energy —eEr balances the force of Coulomb attraction (note that here the 
zero of potential energy is measured at infinite separation). 

Onsager's theory would be expected to be very general and realistic. It should be appli- 
cable, by analogy, to emergent magnetic charges and their transport. In analogy with Eqn. 
[2j the formation of magnetic charges may be represented by the successive equilibria: 

[solvent] = [bound charges] = [free charges] (5) 

With the correspondence e — > Q, E — > B, e — > /xq 1 , we find that the dissociation equilibrium 
into free magnetic monopoles is then characterised by Eqn. [3] with 

6= 8rfTV (6) 

where \x r is the relative permeability of the magnetic 'solvent'. A pictorial representation of 
the physical content of the theory is shown in Fig. 1. 

We define n&, n u as the number of bound and unbound pairs respectively, n Q = rib + n u as 
the total pair concentration and a = n u /n as the degree of dissociation. The dissociation 
constant is given by 

a 2 

K = n -°^. (7) 
1 — a 

The increase in this equilibrium constant with increasing magnetic field (the Wien effect) 
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defines a corresponding change in magnetic moment per unit forward reaction: 

fdlnK\ fdAG°/kT\ A// 
~ V 9B ) ~ W 



(8) 



where AG is the change in Gibbs energy per atom. Using Eqn. [3] we find, for the weak 
field limit, 

A/^=^-, (9) 

which is field-independent. 

Following Onsager we now assume that the reaction rates of the first equilibrium are much 
faster than those of the second equilibrium, in which case all molecules may be considered 
as bound pairs'^. Recalling that a < 1, we see that no ~ ~ N, so the total number 
of defects is approximately constant at a given temperature. Following a small disturbance, 
the relaxation of Aa back to its equilibrium value is determined by charge recombination. 
Onsager showed that the decay is exponential with time constant v = 2{i Q n. where k is the 
conductivity, which is proportional to the equilibrium a. Following Eqn. [9] a fluctuation in 
a causes a proportionate fluctuation in magnetic moment, which decays at the same rate 
v. Thus, measurement of the magnetic moment fluctuation rate as a function of field is 
equivalent to the observation of the magnetic conductivity, and gives direct access to the 
Wein effect: 

v(B) k(B) o(S) I K(B) _ b b 2 

i/(0) k(0) a(0) \ K(0) 2 24 ' { } 

where b is linear in field (eqn. [6]). This equation should be valid provided that both a and 
its change in field Aa are sufficiently small. It may be used to derive Eqn. [T] in the case 
that u r = 1. 
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3 Application to a Real Material 

In spin ice materials like Ho2Ti207 and Dy2Ti207, the magnetic charges or monopoles are 
predicted to be a consequence of the many body nature^ of the dipole-dipole interactions in 
these materials^. In detail, the Ising-like Ho or Dy moments ('spins') are equivalent to proton 
displacement vectors in water ice^. The spins populate the vertices of a lattice of linked 
tetrahedra, pointing either 'in' or 'out' of any particular tetrahedron. The configuration 
'two spins in, two out' corresponds to a water molecule, H2O. At low temperature, magnetic 
interactions of mainly dipolar character ensure a ground state consisting only of such 'two 
in, two out' configurations: the magnetic system is governed by an ice rule and shares with 
water ice the Pauling zero point entropy^. The elementary excitation out of the spin ice state 
is a single spin flip, which breaks the ice rule on neighbouring tetrahedra to create a pair 
of defects, 'three in one out plus one in three out'. The equivalent water ice configuration 
is a closely bound ion pair HaO+OH" 11313 . Castelnovo et a/. m predicted that defects should 
unbind and behave as free magnetic monopoles (in the H-field) that interact via the magnetic 
Coulomb law. This interaction is determined only by fundamental constants except that the 
charge takes the value ±Qtheoreticai = 2/x/a, where fi ~ 10/^^ is the rare earth moment 
and ad = 4.3356A is the distance between tetrahedron centres. The formation of magnetic 
charges in spin ice is thus directly analogous to Eqn. [2J As water ice shows a Wien effect^, 
it seems valid to test for one in spin ice. 

We studied Dy2Ti207 by the technique of transverse field muon spin rotation, with the 
aim of testing the theory of Section [2j and of measuring the 'elementary' charge Q. In /iSR, 
muons implanted into a sample precess around the sum of the local and applied fields, and 
their decay characteristics give information on the time dependence of these fields. In an 
applied transverse field the muon relaxation function has an oscillatory form, resulting from 
the uniform muon precession about the applied field. However this uniform precession can 
be dephased by fluctuating local fields that arise from the sample magnetization M(r, t). 
In the low temperature limit of slowly fluctuating magnetization the dephasing may lead 
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to an exponential decay envelope of the muon relaxation function in which the decay rate 
A depends only on the characteristic rate v of the magnetic fluctuations. By dimensional 
analysis A oc u, so the key property v(B)/v{$) = A(I?)/A(0) can be directly measured (see 



Eqn. 10). In the opposite (high temperature) limit of fast fluctuations, A also depends 
on the width of the field distribution a and the muon gyromagnetic ratio 7. Here, one 
expects A oc (o^) 2 /v (RefJ^), with a independent of temperature, so A is proportional to 
the relaxation timescale r = v~ l . 

The muon relaxation rate \{B) was measured as a function of field and temperature 
after zero field cooling of the sample to 60 mK (see Methods). Selected experimental results 
through the region of interest are shown in Fig. [2] To compare with theory it is useful to 
define an effective magnetic charge by Q c r = 2.1223 ^experimentai^ 2 ^ 3 ( see Eqn. 1), clearly 
demonstrating that only relative changes in A are required to obtain an absolute measurement 
of Qefi- Here we set \i T = 1, which is appropriate to spin ice under the conditions of interest 
(see Methods). 

We would expect Qeff to approximate the true Q only in a finite range of temperature. 
At too low temperature (say b = 2mB ^> 3) the theory breaks down because Aa becomes 
large, while at too high temperature, it breaks down because a becomes large. These limits 
on the range of validity can be estimated experimentally. 

We first consider the temperature dependence of \(B = 2mT), as shown in Fig. [3j The 
high temperature behaviour mirrors the known behaviour of the magnetic relaxation time, 
decreasing with temperature as e constant / T (Ref. I 10 * 19 *). At low temperature A increases with 
temperature, consistent with the expected A oc v oc a. From the graph we can determine 
that the crossover from the presumed unscreened regime to the screened regime with rapidly 
increasing a is at T up per ~ 0.3 K. The inset of Fig. [3] shows the temperature dependence of 
\{B = 2mT, lmT): the apparent irregularity of the measurements below T < T uppei . is in fact 
consistent with the data collapse predicted by theory (see below), and naively demonstrates 
that the current is dependent on B. By measuring b with Eqn. [TOl we estimate that the 



theory starts to break down below Xi ower = 0.07K. Thus 0.3 < T < 0.07 K is the range 
where the theory should be valid (it should be emphasised that this represents a large range 
in the more physically relevant parameter 1/T). 

By fitting the experimental X(B) to extract slope and intercept, we estimate the effective 
charge Qes(T), as described above. In the regime of validity, the value of Q e g saturates 
at ~ 5/is/A, even though both slopes and intercepts depend quite strongly on temperature 
(Fig. [2]). The (inverse) temperature dependence of Q e s is illustrated in Fig. 4, where the 
effective charge can be seen to increase rapidly above T upper = 0.3K, and to decrease slightly 
below Ti ower = 0.07 K. This direction of this slight deviation is consistent with the theory. 

4 Discussion and Conclusions 

We now compare these results with the properties, known and predicted, of Dy 2 Ti 2 07. As 
mentioned above, the predicted monopole charge for this material is ^theoretical = 2/i/a = 
4.6 yU^/AtS Fig. [4] illustrates how the experimentally measured charge very closely ap- 
proaches this value in the temperature range Ti ower < T < T uppcr . The actual value of 
2u PP er may be rationalised on the basis of our expectations for spin ice. Thus the monopole 
concentration should increase as c x e _2J//T , with c C 1 meaning that in the above tem- 
perature range it should vary between, at most, 10~ 14 and 10 -4 . This is indeed a typical 
degree of dissociation for a weak electrolyte that obeys Onsager's theory (i.e. screening is 
negligible)™ 31 . 

The success of the method is impressive. For example, regarding Fig. 2, this data may 
be collapsed over a significant range of temperature and field using a single parameter that 
is in excellent agreement with that calculated by a microscopic theory. However, even more 
remarkably, the method produces precisely the same data collapse (modulo the value of the 
charge) over a range of experimental data on a broad variety of liquid and solid electrolytes^. 
Thus using Onsager's theory, we have demonstrated a perfect equivalence between electricity 
and magnetism. It is further a testament to Onsager's genius that his theory is not only 
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correct and universal, but also robust, in that it may be applied to straightforward experi- 
mental measurements of only modest precision 12 . As already mentioned, the robust quality 
of the method stems from the fact that it deals with relative rather than absolute quanti- 
ties. In our case, this means that the muon spin rotation measurements are not subject to 
interpretation: many different analysis protocols would yield essentially the same result. 

Our /iSR measurements are in full agreement with other /iSR studies of Dy 2 Ti 2 07'221 as 
well as bulk relaxation measurements^. Using ultra low field magnetization measurements 
we have independently verified that a spin ice freezing occurs at around T = 0.7 K. Below 
this temperature the bulk magnetization responds only extremely slowly to an applied field, 
but there are dynamics arising from the magnetic charges. Here we have proved the detailed 
origin of the dynamical response discussed in Ref. ^ 

We have thus established overwhelming evidence that deconfined magnetic charges exist 
in a spin ice material, that they interact by Coulomb's law and are accelerated by an applied 
field. Our experimental results are consistent with those discussed in Refs.tU^, but more 
directly probe the quantities of interest in the regime of temperature and field where a dilute 
gas of magnetic monopoles can be interrogated. We have also measured relative changes 
in the magnetic conductivity (Fig. 2), and proved that monopole currents exist. This 
demonstration raises an interesting question. As discussed by Jaubert and Holdsworth^, 
the dynamics of spin ice monopoles are to some extent constrained by their Dirac strings, 
and this almost certainly precludes a direct current. It is therefore notable that diffusion 
currents do exist, as this would suggest that macroscopic alternating currents may also be 
achievable in spin ice. 

In the broader context, it has been suggested that non Ohmic transport characteristics of 
the Wien effect make it important in biological systems^. Therefore one might imagine that 
there could be novel technological applications for the magnetic equivalent. This is clearly a 
long way off as the temperature scales are low, but it would be interesting to establish proofs 
of principle in this regard, and search for other systems (perhaps micromagnetic arrays^ 
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that show similar effects. In this context, we note the universality of the current method, 
and recommend it as a diagnostic tool in the study of magnetic charge transport. 
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5 Methods 

High quality single crystals of Dy2Ti20"7 were grown using the floating zone method, and 
aligned using an X-ray Laue camera. In the experiment, the crystals were aligned with the 
initial spin polarisation of the muons along the [110] axis, with the magnetic field applied 
at right angles to this axis. In fiSR experiments, one measures the asymmetry of the muon 
beta decay as a function of time A(t), which is proportional to the time evolution of the 
muon spin polarization. A(t) depends on the distribution of internal magnetic fields and 
their temporal fluctuations. 

Our experiment was performed on the MuSR spectrometer at the ISIS facility in the UK. 
All data was taken upon warming after zero field cooling (< 10/xT) from 4K to 60 mK, the 
temperature dependence being repeated twice as a control. Transverse field spectra were 
fitted out to long times using A(t) = Aexp(—[Xt}) cos(2Tivt). The frequency of oscillations 
(v) can be expressed by v = / -f fl \B\/2ir, where B is the average magnitude of the local field at 
the muon site and 7 M is the muon gyromagnetic ratio. A constant, quick gaussian relaxation 
(A = 3/iS 1-1 ) was required for all the data, as expected for spin ice materials^. The fitting 
procedure allowed the slow oscillatory relaxation A (B), and the applied field to vary, meaning 
we are directly reporting the relative change in A. Fig. [5] shows the time evolution of data 
taken at lOOmK in fields of both 1 and 2 mT along with the fits to the data; the envelope 
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of the relaxation is also shown. Background measurements have also been performed on a 
background silver plate to test for the field inhomogeneity in the coils: here the observed 
relaxation was an order of magnitude slower than in the Dy 2 Ti 2 07 sample. 

When zero field cooled below 0.7 K, the susceptibility of Dy 2 Ti 2 07 is sufficiently small 
to approximate /i r — 1 and to ignore demagnetizing effects. 
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Figure 1: Illustration of the magnetic Wien effect, and the detection of magnetic charge by 
implanted muons. (a) In zero field magnetic charges occur mainly in tightly bound pairs, 
but some thermally dissociate to give a magnetic moment (green arrow) that fluctuates due 
to Brownian motion, (b) The competition of Coulomb potential energy (V ~ 1/r) and 
field energy —QBr creates a finite energy barrier to dissociation: the stronger the field, the 
lower the barrier and more dissociation occurs, (c) On application of a transverse field, many 
pairs dissociate as individual charges are accelerated by the field, (d) After applying the field 
the pairs remain separated while more bound pairs form to restore equilibrium. Magnetic 
moment fluctuations due to free charges (green arrows) produce local fields that are detected 
by an implanted muon (/U + ). 
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Figure 2: A quantity that is proportional to magnetic charge conductivity, illustrating a 
non-Ohmic increase in conductivity with field. The measured muon relaxation rate X(B) is 
plotted as a function of field at selected temperatures. Linear fits to the data (lines) are 
shown, from which it is possible to obtain the intercept, slope and consequently the effective 
charge Q eS {T). 
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Figure 3: The temperature dependence of the muon relaxation rate A measured at an 
applied magnetic field of 2 mT. The high temperature regime follows the expected activated 
behaviour as described in the text (fit is shown as a solid line). At low temperature A is 
proportional to the monopole concentration. Its rapid increase above T upper = 0.37^ marks 
a crossover from the regime of weak screening to strong screening of the charges. The inset 
shows the main figure on a logarithmic scale along with data taken in lmT. The fact that 
the apparently irregular behaviour at low temperature is encompassed by the theory helps 
demonstrate that we are indeed measuring magnetic charge transport. 
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Figure 4: Experimentally measured 'elementary' magnetic charge Q e $ in Dy2Ti207. On- 
sager's theory is valid in the regime Ti ower < T < T U p per where the magnetic charges are un- 
screened. The horizontal green line marks the theoretical prediction of Castelnovo, Moessner 
and Sondhi^. 
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Time (\iS) 

Figure 5: The muon relaxation as a function of time at 100 mK in an applied field of 
lmT (solid circles) and 2 mT (open circles), along with the respective fits to the data. The 
envelope of each function is also shown, clearly showing a difference in relaxation at long 
time. 
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